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ABSTRACT 

We report upper limits to the radio and X-ray emission from the newly dis¬ 
covered ultracool dwarf binary WISE J104915.57—531906.1 (Luhman 16AB). As 
the nearest ultracool dwarf binary (2 pc), its proximity offers a hefty advan¬ 
tage to studying plasma processes in ultracool dwarfs which are more similar 
in gross properties (radius, mass, temperature) to the solar system giant plan¬ 
ets than stars. The radio and X-ray emission upper limits from the Australia 
Telescope Compact Array (ATCA) and Chandra observations, each spanning 
multiple rotation periods, provide the deepest fractional radio and X-ray lumi¬ 
nosities to date on an ultracool dwarf, with log jL\^o\)[Hz~^] < —18.1 (5.5 

GHz), \og{L,^y/Uo\)[Hz-^] < -17.9 (9 GHz), and log (L,/Lboi) < -5.7. While 
the radio upper limits alone do not allow for a constraint on the magnetic held 
strength, we limit the size of any coherently emitting region in our line of sight to 
less than 0.2% of the radius of one of the brown dwarfs. Any source of incoherent 
emission must span less than about 20% of the brown dwarf radius, assuming 
magnetic held strengths of a few tens to a few hundred Gauss. The fast rota¬ 
tion and large amplitude photometric variability exhibited by the T dwarf in 
the Luhman 16AB system are not accompanied by enhanced nonthermal radio 
emission, nor enhanced heating to coronal temperatures, as observed on some 
higher mass ultracool dwarfs, conhrming the expected decoupling of matter and 
magnetic held in cool neutral atmospheres. 

Subject headings: brown dwarfs — stars: activity — stars: coronae 


Introduction 


Recently, the discovery of a brown dwarf binary only 2 parsecs from the Sun was an- 
noun ced, making it the third closest system after the Alph a-Gentauri system and Barnard’s 
star flT;uhmanl 120131) . With L7.5 and TO.5 spectral types flBurgasser et al.ll2013b|) the Luh¬ 
man 16AB system (also known as WISE J104915.57—531906.1) has quickly become a bench¬ 
mark for the study of ultracool atmospheres. A unique feature of the Luhman 16AB sys¬ 
tem is the la rge amplitude pho tometric variability (11% in an i z hlter with a period of 
4.87±0.01hr; iGillon et al.l 1201311 of the TO.5 component, with rapid evolution of the global 


•^Center for Exoplanets and Habitable Worlds, The Pennsylvania State University, University Park, PA 
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weather patterns on timescales of abont a day fICrossfield et al.ll2014j) . The two components 
of Lnhman 16AB are separated by 1.5”, or 3 AU flTjnhma.nl l2013h . so not able to influence 
each other via magnetic interactions. 

In principle, photospheric features causing photometric variability could be magnetic- 
or cloud-related; while there have been sporadic measurements of magnetic activity in mid-L 
and later spectral type dwarfs, no clear trends have emerged which connect photometrically 
variable and magnetically active ultracool dwarfs. Since the atmospheres of ult racool dwarfs 


are increasingly neutral, th ey are less likely to support cool magnetic spots (iGelino et ah 


2 OO 2 I: iMohantv et al.l l2002l) than their earlier type stellar counterparts, and the observed 


variability in late L- and T-dwarfs has been attributed to the pres ence of patchy clouds 
flAckerman fc Marled 1200 ll: iBurgasser et al.ll2002l: iMarlev et al.ll2010l) . Despite this, there is 
evidence that in at least some cases, photometric variability of ultracool dw arfs is linked to 
magnetic activity (IClarke et al.ll2003l: iLane et al.ll2007l: iHarding et a.l.ll2013[l 


Magnetic activity signatures in ultracool dwarfs are rare: to date, only 5 L dwarfs and 
one T dwarf have been detected in the radio band, and only 1 L dwarf and no T dwarf has 
been detected in the X-ray band, with a wide range of behaviors displayed among the small 
number of detections. Studies have shown that for late-M to eyly-L dwarfs faster rotation 
results in an increased radio detection fraction flMcLean et al.l l2012h while X-r ay emission 
seems to be suppressed leading to a sort of super-saturation flBerger et al.ll2010l) . 


Dynamo models explain the generation of magnetic helds in ultracool dw arfs by extrapo¬ 
lating convection-driven geodynamo models with strong density stratihcation (IChristensen et al 
20091 1: part of the convected energy flux is converted to magnetic energy to balance ohmic 
diffusion. Such scaling laws predict quite strong magnetic helds, of order 1 kG for a 1 GY 
old, 0.05 Mq brown dwarf with Teg =1500 K and average density of 90,000 kg m“^. These 
scalings do not, however explain why only a handful of ultracool dwarfs of spectral type L 
and T have been detected through radio observations (implying held strengths compatible 
with these extrapolations) while other objects have considerably lower upper limits. Other 
parameters must govern the generation of radio emission and/or held strength. 

Due to the proximity of the Lnhman 16AB system, its magnetic activity can be probed 
with unprecedented sensitivity. An absence of activity signatures would support the pre¬ 
vailing view that large amplitude photometric variability of early T-dwarfs is not connected 
to magnetism, but rather is a consequence of patchy cloud coverage. We report on two 
epochs of radio observations of Lnhman 16AB with the Australia Telescope Gompact Arr^ 
in March and May 2013 and on a Chandra X-ray pointing carried out in November 2013 q 


^We are also aware of observations of Luhman 16AB with the South African KAT-7 array at a wavelength 




















































4 


These observations provide the most sensitive constraints to date on the radio and X-ray 
emission from ultracool dwarfs. 


2. Observations 
2.1. Radio Observations 


Table 1: Upper Limits from Radio Observations 


Frequency 

GHz 

Epoch 

YYYY-MM-DD 

tint 

hrs 

a 

/iJy 

beamsize 

”x” 

pos. angle 
° E of N 

5.5 

2013-03-09 

9 

6.5 

2.3x1.2 

7 

5? 

2013-05-02 

6.5 

8.5 

2.7x1.3 

-28 

5? 

both epochs 

15.5 

5 

2.1x1.3 

-5 

9 

2013-03-09 

9 

7.5 

1.4X0.8 

7 

5? 

2013-05-02 

6.5 

15 

2.0X0.8 

-19 

5? 

both epochs 

15.5 

6.8 

1.4X0.8 
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Luhman 16AB was observed twice with the ATCA: in the 6A conhguration (baselines of 
0.337-5.94km) on 09 March 2013 and again with the 6C conhguration (baselines of 0.153- 
6.0 km) on 02 May 2013 (UT). Continuum mode observations were taken on both dates in 
dual-sideband m ode simultaneously at 5.5 GHz and 9.0 GHz. The Compact Array Broadband 
Backend fCABB iWilson et al.ll201ll) was used with 2 GHz bandwidth per observing frequency 
and 2048 channels each 1 MHz wide. The gain calibrator QSOB1036—52 was used for both 
epochs, with primary hux calibrator QSO B1934—638; Luhman 16AB was tracked in 10- 
minute intervals for both epochs. The hux calibrator was used also as the bandpass calibrator 
during the March observations, but RFI at early times during the May observations prevented 
it from also being used as a bandpass calibrator; instead a single scan of QSOB1036—52 
taken at high e levation was u sed for bandpass calibration. All data were reduced using the 
AIPS package fjCreisen 2003) and best practices for wide-band data reduction. Table [U lists 
the beam sizes for each epoch and frequency band. 


of 20 cm, which took place over 7 hours on April 8, 2013 (R. Fender, private communication). The rms 
image noise is 4 mJy and there is no source near the expected position of the brown dwarf binary. Because 
of the factor of ^1000 disparity in upper limit of the decimeter wavelength observations compared with the 
centimeter wavelength upper limits, we concentrate on the ATCA results in the following discussion. 













- 5 - 


-53 18 45h 


19 00 


15 


30 




_L 




_ ^^^ 

15 14 13 

Right Ascension (J2000) 
Levs = 5.000E-06 * (-3, 3, 4, 5, 6, 22, 38, 55) 


10 49 17 


16 


_L 

12 


Fig. 1.— Plot of the 5.5 GHz image around the expected position of LuhmanlGAB from 
combining both epochs of observations. The box at the center has a length of 3.2”, corre¬ 
sponding to a combination of position uncertainties; see text for details. The beam is shown 
in the lower left of the image. 
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Data at 5.5 and 9 GHz for each epoch were imaged separately and after being combined 
into a single data set per band; both Stokes I and V were searched for Lnhman 16AB. We 
pr opagated t he W ISE position at epoch MJD 55380.018731 by the proper motions given 
in iLuhmaru (120131) to get the expected position coordinates (J2000 10 49 15.99 —53 19 
04.9 for March and 10 49 16.01 —53 19 04.9 for May). There is only a 0.25” difference in 
the position between the two radio epochs. The error range for the expected position of 
Lnhman 16AB is comprised of ±0.75” from the binary separation, an npper limit of ±0.5” 
from parallactic motions (likely smaller dne t o the two m onth separation of the epochs), 
±0.07” from position uncertainty (t aken froniTnhrnanll20131) . and ±0.02” from propagating 
proper motion uncertainties stated in Tankman (12013 ). No source was found near the expected 
position for the brown dwarf pair in either epoch. Figure [1] shows the 5.5 GHz radio sky 
around this expected position, with a box of ±1.6” encompassing the maximum of all the 
errors stated above. The nearest statistically detected source is 12.5” away, with a flux 
density at 5.5 GHz of 30 /iJy. No bursty emission in Stokes I or V is evident in either of the 
bands in either epoc h, in light curves with bin sizes of 60, 300, and 600 seconds (see e.g.. 


Osten fc WolkI 120091) . Details of the observations (and sensitivities derived from individual 


and combined epochs) are listed in Table [H A la upper limit of 5 (6.8) /rJy/beam at an 
observing frequency of 5.5 (9.0) GHz, and at a distance of 2pc, translates into a 3cr radio 
luminosity upper limit Ljy of 7.2x10^° (9.8x10^°) erg s“^ Hz“^ for Lnhman 16AB. 


2.2. X-ray observations 


Ghandra observed Lnhman 16AB on 10 Nov 2013 for 50 ks (ObsID 15705) using the 
AGIS-S3 detector. The data analysis was performed with the GIAO software packag^ version 
4.6. The analysis started with the level 1 events hie provided by the Chandra X-ray Genter 
(GXG). In order to optimize the spatial resolution, pixel randomization was removed. The 
events hie was hltered on event grades (retaining the standard grades 0, 2, 3, 4, and 6), and 
the standard good time interval hie was used. 


We determined the position o f Lnhman 16AB at the epoch of the Chandra observations 
using the proper motion given bv iTaihmanl (120131) . The predicted coordinates are J2000 10 
49 14.65 -53 19 05.2, and there is no evidence a of source at this location (see Figure |2]). As 
discussed in Section 2.1, the uncertainty in the expected position is a maximum of 1.34”. For 
an on-axis source, 90% of the encircled energy lies within a radius of 2” (Ghandra Proposer’s 


^CIAO is made available by the CXC and can be downloaded from 
http: / / CXC . harvard. edu/ciao/download / 
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Fig. 2.— Figure showing the region of the X-ray sky around the expected position of Luh- 
man 16AB at the time of the Chandra observation. Red circle indicates the 90% encircled 
energy radius of 2” centered at the expected position of Luhman 16AB; as discussed in Sec¬ 
tion 2.2, this encompasses the maximum deviation from the expected position by propagating 
errors in position, proper motion, parallax, and binary separation. 






Observatory Guide), and we used this as our guide for establishing the spatial region in 
which to probe for any X-ray emission. The closest detected X-ray source is separated by 
57" from this position, with a total of 17 counts. 


Calculation of an upper limit proceeded with estimation of the background rate at the 
position of Luhman 16AB. 90% of the encircled energy lies within 2" of the central pixel at 
an energy of 1.49 keV (Chandra Proposer’s Observatory Guide). Using an annulus extending 
from 2" — 10" around the position of Luhman 16AB, we calculate a mean background rate 
within 2" of the target, or 0.37 counts (0.2-2 keV) in 48.35 ks. We calculated the quantile 
distribution for a Poisson dist r ibutio n with this intensity in the R statistical computing 
software package IR Core TeamI 020121) , and hnd an upper limit of 2 counts at a signihcance 
level of P=0.001; th is corresponds to a conhdence level of 99.9%, equivalent to a Gaussian 
sigma level of 3.09 (Gehrels 19861) . For the on-source exposure time of 48.35ks, the upper 


limit count rate is then 4.1 ■ 10“® cts/s. We calculated the count-to-flux conversion factor for 
a one-temperature thermal plasma (APEC model) with PIMM^ and we verified that it is 
insensitive for reasonable assumptions on the plasma temperature considering the negligible 
absorption expected for the 2 pc distance of Luhman 16AB. We thus constrain the X-ray 
luminosity in the 0.2-2 keV band to log Lx [erg/s] < 23.0. 


3. Discussion 

3.1. Magnetic Activity Constraints for Luhman 1 GAB 

The upper limits of Luhman 16AB presented in this work are the strongest constraints 
obtained so far for the radio and X-ray luminosity of any ultracool dwarf. We compute the ra¬ 
dio and X-ray activit y indices, log (L r ,%Lh ni) and log (Lx/Lboi), making use of the bolometric 
luminosities given by iFahertv et al.l (1201411 for both components of the binary. We evaluate 
the activity indices separately for the L7.5 and the TO.5 component, assuming that only one of 
the binary components is possibly magnetically active. However, the bolometric luminosities 
of Luhman 16A and 16B are almost the same and the error we make by using their average 
is likely smaller than the sum of all other uncertainties. We hnd log (Li-^iz/Lboi) < —18.1 (5.5 
GHz), log (Lr^jy/Lboi) < —17.9 (9 GHz), and log (Lx/Lboi) < —5.7. Figure [3] puts these upper 
limits in the context of detections and upper limits for other ultracool dwarfs. 

The upper limits at 5.5 and 9 GHz for Luhman 16AB from ATGA are a factor of 15 
and 10, respectively, lower than the most sensitive upper limit of any other ultracool dwarf. 


•^http: / / cxc.harvard.edu/toolkit / pimms.jsp 
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Spectral Type 


Fig. 3.— (Top) Radio luminosity versus spectral type for ultracool dwarfs with spec¬ 
tral types of M6 and later. The bulk of the data refers to the 8 GHz band and is shown 


McLean et ahl ( 

2012 

), 

Williams et ah 

(2012 

) and 

Williams et ah 

(2013 

)• r 


hlled circles correspond to detections. Dotted vertical lines connect measurements of the 
same object; a star symbol connects flare measurements with measurements or limits on 
quiescence of the same object at the same wavelength. The radio upper limits for the Luh- 
man 16AB system are assigned to either the L7.5 or TO.5 component of the system and 
connected by horizontal lines. (Bottom) X-ray luminosity versus spectral type fo r stella r 
objects wit h spec tr al types of M6 and later. Data are tak e n from IWilliams et al.l (12014] ). 
Cook et ah 1 2014 1. Audard et ah (2005), and Stelzer et ah ( 2012 1. Symbols are as in the 


top panel. 
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The X-ray upper limit is about a factor of two deeper than previous sensitive upper limits, 
e Ind Bab (T1-I-T6) is the next c l osest brown dwarf with sensitive limits on X-ray and 
radio emission, from lAudard et ahl (1200511 . At a distance of 3.6 pc, it is only a factor of 1.8 
further away than Luhmanl6AB, conveying about a factor of three difference in luminosity 
sensitivity, and measurements were made for both systems with the same radio and X-ray 
facilities: ATCA and the Chandra X-ray Observatory. The disparity in radio upper limits 
can largely be attributed to the increase in bandwi dth available with th e CABB on ATCA 
now (2 GHz) compared with what was available for lAudard et al.l (1200511 ’s observation (128 
MHz). 


3.2. Interpretation 


We have presented the most sensitive upper limits on X-ray and radio emission for 
ultracool dwarfs to date. The T dwarf component has a measured rotation period of 4.87 h 
(iGillon et al.ll201311 . Our X-ray limit for Luhman 16AB confirms for late L and T dwarfs the 
previous evidence gained from M/L dwarfs for a sharp drop of X-ray activity levels despite 
fast rotation. This absence of X-ray activity is most likely associated with the high electrical 
resistivities in such cool atmospheres which prevent the coupling of matt er and magnetic 


held w hich is necessary to develop magnetic activity (iMohantv et al.ll2002l) . iFleming et ah 


(119951 ) detected stellar coronal heating efficiencies (as measured by Lx/Lboi) down to levels 
approximately a factor of hve lower than our upper limit. Our upper limit is also consistent 
with the activity levels of the Sun at the highest points of its activity cycle, which reaches a 
maximum of logLx/L^oi = —5.9 in the 0.2-2.4 keV band (jPeres et al.l 1200011 . 


As Figure [3] shows, there is a marked drop-off in the number of radio detections for 
objects later than mid-L, with a range of radio luminosities observed at a hxed spectral type. 
Detections of radio emission in ultracool dwarfs are often used to argue for the existence of 
strong magnetic helds. However, the inverse is not true: the lack of radio detection does 
not allow for a determinati on of magneti c held strengths in either of these objects, contrary 
to the statements made in iBergeii (120061 1 for upper limits on radio emission to a sample of 
ultracool dwarfs. Important conclusions regarding the physical extent of any emission can 
be drawn from the radio hux density upper limits in examination of the conditions under 
which these mechanisms operate. 

The interpretation of the variable radio emission in ultracool dwarfs has centered around 
the action of an elect ron-cyclotron maser operating in a region of high magnetic held strength 
(iNichols et ahl l2012l) . In this scenario, the observing frequency is tied to the electron- 
cyclotron frequency in the emitting region and is related to Uc = 2.8 x 10® B (MHz) by 
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t'obs = sUc for harmonic number s equal to 1 or 2, implying kG fields in the radio-emitting 
region detected at cm waveleng ths, and consistent wit h extrapolations from convection- 
driven geodynamo scaling laws fichristensen et ah 2009 ). The intensity of radio emission 
expected from a coherent process such as this is not predictable based solely on the number 
of emitting particles or magnetic field strength. However, the high brightness temperatures 
requi red for coherent emission (usually taken to be Tf, >10^^K: iKellermann fc Paulinv-Toth 


19691 1 coupled with the radio flux density upper limit and observing frequency sets a stringent 


upper limit on the size scale of any radio-emitting region. Rewriting the standard equation 
flDulkl 1198511 for parameters applicable to the current case, an d taking the dwyf ra dius to 
be approximately 1 Jupiter radius, in line with measurements fjSorahana et al.ll2013h . leads 
to the following constraints: 


xr, = 1316- 


d 


pc 






1/2 


( 1 ) 


where x is the size of any radio-emitting region in units of Rj {Rj is one Jupiter radius 
(=7.1x10® cm)), dpc is the distance to the dwarf in pc, uomz is the observing frequency in 
GHz, S^jy is the flux density in /xJy, and Tf, is the brightness temperature in K. Evaluated for 
the upper limits at the two frequencies (and assuming =10^^K) gives x <0.002 at 5.5 GHz, 
and X <0.001 at 9 GHz. These are upper limits, as stellar phe nomena have demonstr ated 
the existence of brightness temperatures as high as Tf, r^IO^^K fjOsten fc BastianI 1200811 . 


The upper limit on size holds if the conditions are right to produce coherent emission. 
Growth rates of the cyclotron maser instability are maximized in relatively rarehed, magne¬ 
tized plasmas where the dim ensionless ratio of the plasma frequency to the ele ctron-cyclotron 


freque ncy is less than a few flLee et al.ll2013fl . The atmosphere calculations of iMohantv et ah 


(120021 ) showed that in the lower atmosphere, the total density for a dusty atmosphere model 


in a cool dwarf with Teg near 1500 K will be about 10“® g cm“^, with an ionization fraction 
of about 10“^^. This would suggest an electron density of approximately 10^ cm“^. The kG 
field strengths for these objects derived from scaling laws, combined with these parameters, 
indicate that the conditions for the instability may exist, but the cyclo t ron m aser mechanism 
could be inoperable for reasons still to be determined. iMutel et al.l (120071) found a strong 
dependence of the growth rates of the cyclotron maser instability on the opening angle of 
the loss cone distribution of electrons that could power the instability. Beaming effects may 
also explain the lack of detections, if there is a misalignment between the opening angle of 
the emission and the line of sight. 

Another possibility that has been put forward to explain the quiescent radio emission 
from ultracool dwarfs is gyrosyn chrotron ern ission, in analogy with the magnetic activity 
seen in higher mass dwarf stars (iGiidell 120021 ). For this incoherent process, the strength of 
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the emission depends not only on the magnetic held strength in the radio-emitting source, 
but also on the index of the distribution of accelerated particles with energy (5), and the 
size of the emitting region. Figure 0] displays the values of 6, B, and size of the emitting 
source that are compatible with the observed upper limit on hux dens ity at the tw o radio 
frequencies, given the limit of applicability of the analytic expressions in lDulkl fjl985l) . While 
the constraints on size are not as stringent as for the case of a coherent emission, they do 
rule out a global gyrosynchrotron-emitting magnetosphere around one of the dwarfs in the 
Luhman 16AB system, as this would lead to detectable levels of gyrosynchrotron emission. 

Our discussion has concentrated mainly on steady-levels of emission. The observation 
of radio bursts in a relatively small sample of all radio-observed UCDs combined with the 
small rotation period of the known radio bursters has given rise to the discussion of selection 
effects, e.g., the typical length of the radio observations (few ho urs) may not have covered 
the - generally unknown - full rotational cycle of many UCDs (jStelzer et ahl 120121) . This 
bias can be ruled out here. Multiple rotation periods of the TO.5 dwarf were covered with 
our radio data, so the observations were sensitive to bursts occurring at particular rotation 
phases. Yet, the viewing geometry and/or topology of the magnetic field may prevent the 
detection of such bursts — if present — on Luhman 16AB. Limits presented here are likely 
unachievable for other ultracool dwarfs in the near future, only to be exceeded possibly by 
measurements from the Athena mission (for X-rays) and the Square Kilometer Array or Next 
Generation Very Large Array future radio telescopes. 
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